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Abstract 

Generalized parton distributions at large x are studied in perturbative QCD approach. As x — > 1 
and at finite t, there is no t dependence for the GPDs which means that the active quark is at the 
center of the transverse space. We also obtain the power behavior: H^{x,^,t) ~ (1 — x) 2 /(l — £ 2 ) 
for pion; H q (x, £, t) ~ (1 — x) 3 / (1 - £ 2 ) 2 and E q (x, £, t) ~ (1 — x) 5 /(l - £ 2 ) 3 /(£) for nucleon, where 
f{€) represents the additional dependence on £. 
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In recentyears, there has been considerable interest in generalized parton distributions 
(GPDs) 0,12,3, 

which were introduced originally to understand the quark and gluon contri- 
butions to the proton spin m . They are also related to the quantum phase space distributions 
of partons in the hadrons [4|. The theoretical framework of the GPDs and their implications 
about the deeply virtual Compton scattering, deeply virtual meson production, and the 
doubly-virtual Compton scattering have been well established 

H BBS SHIS- Apart 

from the renormalization scale, the GPDs depend on the momentum transfer t, the light-cone 
momentum fraction x, and the skewness parameter £ which measures the momentum trans- 
fer along the light-cone direction. In phenomenology, the GPDs are parameterized through 
the double-distributions 12|| and fit to the experimental data 0, H, El Ei EH ■ However, these 
parameterizations have too much freedom, and we still have a long way to go for a complete 
understanding of the GPDs. In this context, any theoretical result on the behavior of GPDs 
will provide important information. For example, the polynomality condition 5], and the 
positivity constraints E3 have already played significant roles in the parametrizations of 
GPDs. The light-cone framework provides useful guidelines for calculating the GPDs once 
the wave functions are known [14J. More recently, the GPDs at large t have been explored 
[Tif . yielding important constraints as well. 

In this paper, we study the GPDs in the kinematic limit of x — > 1. For the forward 
parton distribution, a power behavior at large x was predicted based on the power counting 
rules, for example, (1 — x) 2 for pion, and (1 — x) 3 for nucleon 0, BElI, m H This 



power behavior comes from the fact that the hard gluon exchanges dominate the structure 
functions at x — > 1, and is calculable in perturbative QCD j^Elf. In this paper, we will 
follow these ideas to analyze the dependence of GPDs on the three variables x, £ and t in 
the limit of x — > 1. We use the QCD factorization approach, and express the GPDs in terms 
of the distribution amplitudes of hadrons. In the limit of x — > 1, the power behavior of the 
GPDs does not depend on a particular input of the distribution amplitudes, and therefore 
can be predicted model-independently. More importantly, the £ and t dependences can also 
be calculated. For example, we find that there is no t-dependence at x — > 1, which agrees 
with the previous intuitions 24] . 



We take (1 X^j clS db small parameter, and expand the GPDs in terms of (1 — x). In 
the process, we assume the variables £ and t finite. Finite £ means £ < x and restricts 
our analysis valid in the DGLAP region for the GPDs. The relevant Feynman diagrams are 
shown in Fig. 1 for pion, and in Fig. 2 for nucleon for a typical contribution. The variables P 
and P' are the initial and final state hadron momenta, respectively, and t = A 2 = (P — P') 2 . 
We further introduce two vectors P and n: P = (P + P')/2, n 2 = 0, and n ■ P = 1. The 
skewness parameter is defined as £ = — n ■ (P' — P)/2. The initial and final light-cone 
momenta of the quarks are then {x + £) and {x — £), respectively. In the following, we will 
neglect the masses of the hadrons, and then t = — A 2 L /(1 — £ 2 ) where Aj_ is the transverse 
part of the momentum transfer A. 

As shown in Fig. 1, the intermediate state has momentum k which will be integrated out. 
To avoid an infrared divergence, we keep k± much larger than Aqcd- The offshellness of 
the quark and gluon propagators are on the order of k\/{\ — x). So, we have the following 
hierarchy of scales in the limit of x — > 1: k\j (1 — x) ^> k\ Aq CD , and k\j (1 — x) (— t) 
as well. These relations will be used to get the leading-order results, and any higher power 
in (1 — x) will be neglected. 
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FIG. 1: Leading order contributions to the generalized parton distribution H q (x,^,t) for pion at 
large x. The crosses represent the intermediated states, and the double lines for the eikonal term 
from the gauge link. 



The GPD H(x,^,t) for pion is denned as 
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H q {x^t) =- I —e 
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where L represents the light-cone gauge link. We work in Feynman gauge, and the leading- 
order diagrams were shown in Fig. 1, where the double lines represent the eikonal contribu- 
tions from the gauge link, and the cross indicates the intermediate state on mass shell. The 
initial and final states are replaced by the light-cone Fock component of hadrons with the 
minimal number of partons. After integrating over the internal transverse momentum l±, 
the light-cone wave function leads to the distribution amplitude, 4>(x) = / ^U^ip(x,l±). 
The calculation of the diagrams in Fig. 1 is straightforward, and the result is 



d 2 k\ 



(2vr) 3 (2k ■ P)(2k ■ P')' 
The integral X depends on the distribution amplitudes of the initial and final states, 



(1) 



T(x, £) = (47ra s CV) 2 J dx\dy 
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1 + 



where Cp = 4/3 and x\ = 1 — x±. The integral becomes a constant in the limit of x — > 1. In 
Eq. the denominator factors 2(k ■ P) and 2(k ■ P') in come from the gluon propagators in 
the diagrams. They depend on the momentum transfer t in general. However, if expanded 
at small (1 — x), they become 



2k P 
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1 — x 
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1 — X 



2k- P* fc»(l-OL 4(l-0 2 A;i 



1 + 
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4(1+0^1 

-e)t 



X 



(2) 



which implys that there is no t-dependence in the leading order, and any dependence will 
be suppressed by a factor of (1 — x) 2 . Since the propagators are the only source of the 
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t-dependence of the GPD in Eq. (JI}, we conclude that at x — > 1 the GPD H(x,^,t) for 
pion has no t-dependence, and any t-dependence must be suppressed by a factor of (1 — x) 2 . 
These conclusions agree with the analysis in the impact parameter dependent picture of the 
GPDs at £ = [24|, while our results are valid for any finite value of £. 

Collecting the above results, we have the GPD H(x,^,t) for the pion in the limit of 
x — > 1, 

EHx f t) - {1 ~ X)2 f I j ( 3 ) 

There is an infrared divergence where the transverse momentum k±_ becomes soft. This 
divergence breaks the factorization in principle. However, this does not change the power 
behavior. Nevertheless, we include a regulator A to regulate such divergence |23|. Like the 
forward parton distributions, a power behavior is found here for the GPDs in Eq. (JHJ). If 
we take £ = and t = 0, Eq. (jH} will reproduce the forward parton distribution of pion. 
So, in the limit of x — > 1, the GPD H(x,^,t) for pion can be related to the forward parton 
distribution q*(x), 

1 

1 

We note that the above equality saturates the positivity constraints 0, E3] if we take the 
power behavior for valence quark distribution q n (x) ~ (1 — x) 2 . 

We turn now to the study of GPDs for the nucleon. Since the leading Fock component 
of nucleon has three partons, many more diagrams will contribute. Here we only show a 
particular diagram in Fig. 2. There are two intermediate momenta, ki and k 2 . Similar to the 
above analysis for the pion case, we have a hierarchy of scales: (k\) /(1 — x) ^> (k 2 /) ^> Aq CD 
and {kj_)/(l — x) ^> (— t), where (k 2 ^) represents the typical transverse momentum scale, 
(k 2 /) ~ ~ Again, we are only interested in the leading-order result, and neglect 

any higher-order corrections in — 

The calculations are performed in the helicity bases for the initial and final nucleon states 
(A', A), in which the following off- forward matrix elements are defined: 



H*(x^ 1 t) = —-q*(x) . (4) 



Wva = / = I — e*» / P', X' 



PX 



The helicity non-flip amplitude has contributions from both H and E GPDs, while the 
helicity flip one only has the contribution from E GPDjsf, 



H n = Hu = H q (x,£,t) - ^—^E g (x,U) 

n n - 2^r^'^- n - (5) 



We will show how the diagram in Fig. 2 contribute to these amplitudes. 

The helicity non-flip amplitude for the diagram of Fig. 2 has the following form, 

r d 2 k 1± d 2 k 2± r da f 1 1 V P' ■ h 

~ J (2vr) 3 J a(3(l-x)\2P-(k 1 + k 2 )2P'-(k 1 + k 2 )J P> ■ k 2 U ' [ ) 

where any other factors (such as color factors and coupling constants, etc.) are included in 
the integral I p (x, £). This integral depends on the leading-twist distribution amplitudes of 
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FIG. 2: The typical diagram contributing to the nucleon GPDs at x — * 1. The quark helicity con- 
figuration is {TIT}- The bars indicate the places where we need to consider the internal transverse 
momentum expansion to get helicity flip amplitude. 



the proton [25[, and will become a constant integral at the limit of x — > 1. The longitudinal 
momentum fractions of k\ and k 2 are defined as n ■ ki = a(l — x), n ■ k 2 = (3(1 — x), and 
a + (3 = 1. For the propagators, we make expansions at small (1 — x) as before, for example, 



1 (kl) 



2P ■ (h + k 2 ) 1-x 



+ 



t 



l + 0((l-x) 2 )^- + 
(kl) 



(7) 



which again has no t-dependence at the leading order, and any t dependence is suppressed 
by a factor of (1 — x) 2 . All propagators in Eq. @ have this property, and all other diagrams 
which contribute to TL^ at the leading order have the same dependence on t. The t depen- 
dence of nucleon GPDs is the same as that of the pion: there is no t dependence and any 
dependence is suppressed by a factor of (1 — x) 2 . In addition, every diagram contributes the 
same dependence on £. Adding all of the contributions together, we get 

(l-xf rd 2 k 1± d 2 k 2± [da 
= (73^ J (27r) 3 J A)2 p , (8) 

where the function F is of order 1 at x — > 1. Here we also include a regulator A in F to 
regulate the infrared divergences in the ki± and k 2 ± integrations. If we take £ = and t — 0, 
the above results will reproduce the forward parton distribution at large x. That means we 
can have, 

1 (1-x) 3 

^TT = (1 _ e)2 ^) ~ (l_^2)2 ' ( 9 ) 

at x — > 1. 

Since hard scattering conserves the quark helicity, in order to get the helicity flip am- 
plitude Ti^i we must include non-zero orbital angular momentum either for the initial or 
final states. In other words, we need to consider the li ght- cone Fock components of hadrons 
with at least one unit of orbital angular momentum [26(. The calculation for the helicity 
flip amplitudes are much more complicated than that for the helicity conserving ones. The 
method we are using follows Ref. [27( where the helicity flip Pauli form factor was calculated 
in perturbative QCD. We will sketch the method and summarize the main results, but skip 
the detailed derivations. 

First, we keep the internal transverse momenta l± of the scattering partons in the hard 
partonic scattering amplitudes. Then, we expand the amplitudes at small l±. Since Aj_ 
is the only relevant external transverse momentum, the expansion of the amplitudes will 
be proportional to Aj_ • l± or Aj_ x l±. Integrating these terms over l± with the light- 
cone wave functions, we will get, e.g., / d 2 l±A ± ■ + ill)^ 3 ' 4 ^ ~ (A^ + zA^)<3>( 3 ' 4 ), and 
/ d 2 l ± A ± x + Wl)^ 3 '^ ~ — i(A.j_ + zA^)$( 3 ' 4 ), where ?/>( 3 > 4 ) are the light-cone wave 

functions for the Fock state with one unit orbital angular momentum j26[, and $( 3 ' 4 ) are the 
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related twist-four distribution amplitudes [25] . Thus, the final results of TCn depend on the 
twist-three and twist-four distribution amplitudes of the nucleon. 

We must consider the expansions for all propagators and quark wave functions which have 
dependence on As an example, in Fig. 2 we indicate all places where the l± expansion 
should be considered if the initial state has one unit of orbital angular momentum. These 
expansions will give additional power of (1 — x) 2 , leading to the helicity flip amplitudes 
suppressed by (1 — x) 2 . For instance, one of the gluon propagators in the diagram of Fig. 2 
has the following expansion, 



1 1 



(k 2 -x 3 P-l ± ) 2 (k 2 -x 3 py 



x _ 0(1 - x) 2 A ± ■ h 



2± 



(10) 



Extracting the expansion coefficients, and combing with other factors in the amplitude, we 
get the contribution to the helicity flip amplitude Ti^ from this term: ~ (1 — x) 5 /(l — 
£ 2 ) 2 (1 + £) 2 . Adding the similar contribution from the final state expansion, we get TC^ ~ 
(1 - x) 5 /(l - e 2 ) 2 (l/(l + 2 + V(l - 2 ) = (1 - x) 5 (l + £ 2 )/(l - £ 2 ) 4 - All expansions result 
in the same suppression of (1 — x) 2 . However, they do not contribute the same dependence 
on £. For example, the quark wave function expansions lead to TC^ ~ (1 — x) 5 /(l — £ 2 ) 4 . In 
summary, the helicity flip amplitude will have the following result at x — > 1, 

W 4t ~(AI + iAi)^£^/(0. (11) 

Here /(£) represents an additional dependence on £, which will depend on the input of the 
twist-three and twist-four distribution amplitudes of the nucleon. From this, we deduce the 
behavior of GPD E q (x,£,t) as, 

E q (x,U)~^^m- (12) 

Comparing with Eq. 0, we can neglect the E q contribution to the helicity non-flip ampli- 
tude, and then we have H q (x, £, t) ~ (1 — x) 3 /(l — £ 2 ) 2 . So, in the limit of x — > 1, H q (x, £, t) 
can be related to the forward quark distribution q(x), 

H q (x,£,t)= {1 ^ e)2 q(x) . (13) 

Again this relation saturates the positivity constraint 0, ^| for nucleon GPDs if the forward 
quark distribution takes the power behavior at large x: q(x) ~ (1 — x) 3 . 

Before concluding, a few cautionary comments are in the order. First, we omit the scale 
dependence of the GPDs. The scale dependence at large x is not just the simple DGLAP 
evolution 21, 28|. In our calculations we implicitly assume Q 2 (l — x) ^> Aq CD . Second, at 
the limit of x — > 1 there exist a"log m (l/ (1 — x)) for m < 2n series terms which need to be 
resummed, leading to a Sudakov form factor suppression [i^, H^, IH, |2t| . Third, the soft 
mechanism might contribute to the GPDs |30j at x — > 1. We did not include such effects in 
our analysis. 

In summary, we have studied generalized parton distributions at x — > 1. We found 
that the pion's GPD H*(x,Z,t) ~ (1 - x) 2 /(l - ^ 2 ), and the nucleon's GPD H q (x,£,t) ~ 
(1 - xf/{l - Z 2 ) 2 and E q (x,£,t) ~ (1 - x) 5 /{l - ^ 2 ) 4 /(0- Th ere is no t dependence, and 
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any dependence is suppressed by a factor of (1 — x) 2 . These results can provide important 
information on the GPDs' parameterizations. 

We thank Xiangdong Ji for pointing out the t dependence of the GPDs at large x and 
many other stimulating discussions. We also thank Andrei Belitsky, Stan Brodsky, and 
Jianwei Qiu for useful comments. This work was supported by the U. S. Department of 
Energy via grants DE-FG02-93ER-40762. 
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